Unconventional superfluidity of fermions in Bose-Fermi mixtures 
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We examine two dimensional mixture of single-component fermions and dipolar bosons. We 
calculate the self-enregies of the fermions in the normal state and the Cooper pair channel by 
including first order vertex correction to derive a modified Eliashberg equation. We predict ap- 
pearance of superfluids with various non-standard pairing symmetries at experimentally feasible 
transition temperatures within the strong-coupling limit of the Eliashberg equation. Excitations in 
these superfluids are anyonic and follow non-Abelian statistics. 



PACS numbers: 

I. INTRODUCTION 

Experimental realization of fermionic superfluidity in 
a quantum degenerate ultracold gas llj started a renewed 
interest in the field. Attraction between fermionic parti- 
cles favours pairing of fermions resulting in superfluidity 
of the system. The paired fermions, known as Cooper 
pairs, can have different kinds of internal symmetries. 
The common ones found in nature have s-wave and d- 
wave internal structure and conserve parity and time 
reversal symmetry. In superfluid "^He, the so-called A 
and Ai phases are characterized by Cooper pairs with 
nonzero magnetic orbital momenta. Also, there was a 
theoretical prediction of higher order /-wave pairing in 
^He 0, [|. Cooper pairs with chiral Px + ipj,-wave in- 
ternal structure are believed to be responsible for the 
observed superfluidity of electrons in Strontium Ruthen- 
ate This kind of pairing breaks the time-reversal 
symmetry. Spinless chiral p-wetve superfluid state has 
formal resemblance to the "Pfaffian" state proposed in 
relation to the fractional quantum Hall state with flUing 
factor 5/2 [1-0]. When conflned in a two dimensional 
geometry, excitations in the chiral p-wave superfluid be- 
come so called non-Abelian anyons. Anyons are particles 
living in a two dimensional plane that under exchange 
do not behave neither as bosons nor as fermions. For 
non-Abelian anyons, exchange of two such particles de- 
pends on the order of the exchange [8l-[l0|. Apart from 
fundamental interest in the existence of such particles, 
non-Abelian anyons find remarkable applications in the 
field of quantum information for quantum memories and 
fault-tolerant quantum computation ^llj . 

Recently it has been shown that quasi-particles in vor- 
tex excitations of chiral two-dimensional jp-wave spinless 
superfluids obey non-Abelian statistics [ll, [11] . Using p- 
wave Feshbach resonances in fermionic ultracold atoms, 
such superfluids can be realized in principle, but this pro- 
cedure is very difficult because of non-elastic loss pro- 
cesses Another proposal for p-wave superfluid with 
dipolar fermions has been recently formulated in Ref [l5| 
with transition temperature in the order of .le/, where 
£/ is the Fermi energy. Bose-Fermi mixtures are another 
candidate for creating superfluidity in fermions via boson 



mediated interactions [16| and have formal resemblance 
to phonon mediated superconductivity in metals. It was 
found, however, when the bosons and single-component 
fermions are completely mixed, increasing the boson- 
fermion interaction strength, or fermionic density may 
induce dynamical instability of the condensate resulting 
in phase separation between the mixture jlGt - tlOi] . Near 
phase separation, inclusion of the dressing of phonons 
is predicted to increase the transition temperature. In 
Ref. [11], the authors found that close to Feshbach reso- 
nances, a Bose condensate of dimers can induce a strong 
pairing in the p-wave channel. 

In the present paper, we discuss another way to gen- 
erate high temperature superfluids in a Bose-Fermi mix- 
ture. We study the property of superfluidity in Bose- 
Fermi mixtures, where bosons are interacting via long- 
ranged dipolar interactions. We show that the transition 
temperature for p-wave superfluidity can become compa- 
rable to the Fermi energy. More importantly, we find that 
other more exotic Cooper pairs with /- and h-wetve inter- 
nal symmetries are possible in certain range of Fermi en- 
ergies without bosons and fermions separating. In addi- 
tion, we study the excitations in chiral states of the odd- 
wave superfluids and point out their non-Abelian anyonic 
nature. Experimentally, an available bosonic species, 
where prominent dipolar interaction can be achieved us- 
ing Feshbach resonance is Cr^^ [H, [23| . Another route 
towards achieving ultracold dipolar gas is to experimen- 
tally realize quantum degenerate heteronuclear molecules 
[11], which have permanent electric moment. Thus, in 
the near future a quantum degenerate mixture of dipolar 
bosons and fermions will be achievable experimentally. 

In Section[II]we review the properties of a dipolar Bose 
condensate in a pancake trap. We discuss specially the 
excitation spectrum of such condensates. Then in Sec- 
tion mil we study the boson- fermion interaction and the 
many-body effect of fermions on dressing the excitation 
spectrum of the condensate. In Section IIVI we discuss 
the interaction between the fermions mediated by bosons 
in different angular momentum channel. By integrating 
out the bosonic mode, we show that the absolute value of 
the interaction between fermions in the p, h, /-wave an- 
gular momentum channel are comparable depending on 
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the Fermi energy. In Section |Vl we go beyond Migdal's 
limit and include first order vertex corrections to study 
fermion self-energy in normal state as well as the mass 
renormalization function in the high temperature limit. 
In doing so, we include the full effect of retardation and 
strong momentum dependance of the bosonic excitation 
spectrum and bosonic propagator. We find that vertex 
correction reduces the mass renormalization function. In 
Section |Vl] we calculate the self-energy in Cooper-pair 
channel for p, /, /i-wave order parameters including the 
vertex correction and cross-interaction for temperature 
T > Tc- By deriving a vertex-corrected strong-coupling 
Eliashberg equation, we solve for transition temperatures 
in different angular momentum channels. In Section [VIH 
we present a brief discussion regarding the possible oc- 
currence of non-Abelian Majorana fermions for broken 
time-reversal p, /, and /i-wave superfluids. We solve the 
Bugoliubov-deGennes equation for the p, /, and h-wave 
superfluids in the limit of large distance to find the Ma- 
jorana bound states. 



the Hamiltonian of the dipolar bosons in the condensed 
phase, Hb — J2k'_i_ ^o{k±)/3j^- Pj^^ , where /3^- and /3^- are 

Bugoliubov operators. The excitation frequency rto{k±) 
is given in the units of trap frequency, 

^lik^eo) = +53d|^V (k^io^^ [k^eo?, (2) 

where £q — ^Jhjrn]pj\,. We define a dimensionless dipo- 
lar interaction strength = STrnibgddnbio / ^fi-^ which 
will be used later, where io is the ground state oscillator 
length. Also the phonon propagator Do{iu!s, k±) in this 
regime is given by 



C.2 + J72(fc^)' 



(3) 



where bosonic Matsubara frequency uj^ = 2s where s is 
an integer. By minimizing the mean field energy of the 
Bose condensate within Thomas-Fermi regime we find 
that 



II. DIPOLAR BOSE-EINSTEIN CONDENSATE 

Our system consists of dipolar bosons mixed with sin- 
gle component fermions, confined in a quasi- two dimen- 
sional geometry by a harmonic potential in the z direc- 
tion with the condition mi,uj^ — 'm-fUj'j — where mi, and 
m/ is the mass of bosons and fermions respectively and 
ojbjUJf is the trapping frequency. First, we assume that 
the bosons are polarized along the z direction. The dipo- 



lar interaction reads Vdd 



[ikl/k^ - 1) in mo- 



mentum space, where is the dipole-dipole interaction 
strength. For atoms gdd — /^o/^m/47''i and for dipolar 
molecules ^dd = Me / 47reo , where /im and /Xe are the mag- 
netic moment of the atoms and the electric dipole mo- 
ment of the molecules, respectively. We assume that the 
z dependence of bosonic density is given by a Thomas- 
Fermi profile. 



nb{x,y,z) = 1^ ( 1 



R,/£o = 2.5g3d 



167rgdd 



1/3 



For 



3g 

47rgdd 



> 1, the excitation spectrum of the conden- 
sate, denoted by n{k±), can be divided into two parts: 
i) ^ k±, phonon spectrum for small momenta and ii) 
~ fc^, free- particle like spectrum for higher momenta 
[26| . For < 1 and gad greater than a critical value, 

^{k±) has a minimum at momentum ko [29| . where ko 
is in intermediate momentum regime. Following Lan- 
dau, the excitations around the minimum are called "ro- 
tons" . With increasing g^^ the excitation energy at k^ 
decreases and eventually vanishes for a critical particle 
density. When the particle density exceeds that critical 
value, the excitation energy becomes imaginary at finite 
momentum and the condensate becomes unstable. The 
use of Thomas-Fermi density profile is justified in this 
region as the chemical potential necessary to reach roton 
instability exceeds hujz- 



where the Thomas- Fermi radius Rz is determined varia- 
tionally. After integrating over z dependence of the den- 
sity profile of bosons, the total interaction takes the form 



5R^ 



V{k±) where 



vCk±) = 



3g 

Svrgdd 
15[2~kl 
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3(l + fcj_)2exp(-2fc_L)-f 3] 



U5 
"J- 



(1) 



k± = k±Rz and g is the contact interaction between 
the bosons. V(fc_L) is repulsive for small momentum and 
can be attractive in the high momentum limit depend- 
ing on the contact interaction. Subsequently, we write 



III. BOSON-FERMION INTERACTION AND 
DRESSED EXCITATION 

In this section, first we disccuss the condensate-fermion 
interaction Hamiltonian and dynamical stability of the 
condensate. We find the the dressed propagator for the 
Bugoliubov quasi-particles in presence of fermionic parti- 
cle hole excitation within second-order perturbation the- 
ory. We then discuss the appearance of roton instability 
in the dressed excitation spectrum of the condensate. 

Kinetic energy for the single component non- 
interacting fermions, is characterized by the Hamiltonian 



H 



f 



where c^- 



and c,- 

k 



fermionic creation and destruction operator. S.{k±) 
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A:^/2m/ — e/ is the dispersion energy of the fermions and 
is the Fermi energy. The density profile of fermions 
along the z direction is approximated by a normalized 
Gaussian with width — -y/ft/m/w/. The fermions are 
interacting with the bosons via short ranged contact in- 
teraction of strength gbf- After integrating over the z 
coordinate, the boson-fermion interaction Hamiltonian 
reads 



where the two dimensional Lindhard function 



fcj_) = 1 — |w| 



3gbfQ 



E 



=1- c 7 b\ bqi , 

k± '-± ?! 



(4) 



and uj is the transfer of energy. As we are interested in 
Cooper instability of the fermions which happens for mo- 
menta close to Fermi momentum, the transfer of energy 
is of the order of w ^ ej. In this limit the Lindhard func- 
tion h{Lo,k±) = 1. For future use, we define an effective 
interaction strength between the bosons and fermions, 



where 6^- , 6^- are bosonic creation and annihilation op- 
erator, and a = exp [-jr^ - ^ (-^ - 2^ erf (^^^ , 

with erf(..) being the error function. Next, we con- 
centrate our attention on the interaction between the 
fermions and the Bugoliubov quasi-particles. The fiuc- 
tuations in fcrmionic density couple to the density fluc- 
tuations present in the Bose condensate. Due to the 
momentum dependence of the bosonic excitations, the 
condensate-fermion interaction becomes function of mo- 
mentum. Including the fluctuations in the Bose conden- 
sate in the x—y plane, the condensate-fermion interaction 
Hamiltonian can be written as 



Ghf 



45ggfjVo 
1287rgdd4 



In order to determine dynamical stability of the uni- 
form condensate, we study the properties of the dressed 
spectrum of Bugoliubov quasi-particles in the parameter 



regime of 



3g 

47Tgdd 



> 1, where the bare spectrum has no ro- 



h,- +V - 



(5) 



ton minimum. From Eq. ([7]), due to the attractive effect 
of Bose- Fermi interaction, we notice that roton minimum 
in the excitation spectrum f2(fc^) can develop for certain 
critical interaction strength g^d and Ghf- Consequently, 
the uniform condensate will become instable as the ex- 
citation spectrum becomes imaginary at a flnite momen- 
tum. But, as long as the roton gap remains positive, the 
uniform state of the bosons will be stable or metastable 
[1^. In Fig. [21 we have plotted the dressed excitation 
spectrum for different values of Bose-Fermi interaction 



where the momentum dependent coupling constant is Gm- With no boson-fermion interaction, as 



3g 

47r£/dd 



given by 7(fcj_) = -W 2nbtb(k±) IQ.o{k±). Due to the inter- 
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FIG. 1: Schematic diagram of phonon propagator includ- 
ing the polarizing effect of the fermions. The thick dash- 
dotted line corresponds to the dressed phonon propagator 
D{iu}s,k±), while the thin dash-dotted line corresponds to 
non-interacting phonon propagator Do{iiJs, k±). The solid 
line corresponds to free fermion propagator. 



action between Bugoliubov quasi-particles and fermions, 
as shown in Fig.[Tl the bosonic excitations are dressed by 
the fermions resulting in dressing of bare condensate ex- 
citation spectrum in Eq. ([2]). This can be derived within 
second order perturbation theory as shown in Fig. [T] 
Subsequently the new phonon propagator is given by 



> 1, 

the excitation spectrum has the usual phonon regime for 
low momenta and free particle regime for higher mo- 
menta. For a critical Gbf, i-e. Gbf ~ -5 in Fig. [21 in 
the region of intermediate momenta, the excitation spec- 




FIG. 2: The dressed excitation spectrum Q{k±) as a func- 
tion of k± For various values of the Bose-Fermi interaction 
Gbf ~ 0(blue line), .51(red line), .57(black line). The fixed 
parameters are ruf = 53, -4§f^ = -6, and gzd ~ 2.5 



D{iu!s, k±) = -- 



»o(fc±) 

n^ikA 



(6) 



where the dressed excitation spectrum fl{k±io) for the 
dressed Bugoliubov quasi-particles is expressed by 



(7) 



trum have roton-maxon character. With further increase 
of Gbf , the roton gap goes to zero and the excitation spec- 
trum becomes imaginary at a finite momentum. Subse- 
quently the bosons undergo a phase transition to a state 
with periodic dencity modulation, which finally collapses. 
This periodic density wave state can be stabilized for re- 
pulsive contact interaction obeying Sg/Sirg^d ~ Ghf-^ 
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[28| . The appearance of roton minimum in dressed exci- 
tation spectrum in this case is entirely due to the many- 
body effect of the fermions on Bugohubov quasi-particles. 
Another point we hke to stress is that, the roton instabil- 
ity is always reached before the phonon instability point- 
ing towards phase separation [l6|. 



IV. EFFECTIVE INTERACTION BETWEEN 
FERMIONS 



Next, we look into the variation of Al(0, 0) as a func- 
tion of 77. For concreteness, we assume a ^^Cr-^'^Cr mix- 
ture. The interaction strengths from Eq. ^ in various 
angular momentum channels have been plotted in Fig. [31 
We find that with changing dimensionality, the strengths 
in different angular momentum channel varies. Addition- 
ally, |Ai (0,0)1- |A3(0,0)| - |A5(0, 0)1, and depending on 
the dimensionality they can be positive or negative. 



In this section we find the interaction between the 
fermions mediated by the Bugoliubov quasi-particles in 
the limit of T = 0. We then obtain the interaction 
strength in different angular momentum channel, which 
varies as a function of dimensionality parameter, defined 
as ?7 = ef/hujf. 

Integrating out the bosonic degree of freedom in 
Eq. ([5]), and inclusion of the effect of dressed excita- 
tion spectrum, results in effective interaction between the 
fermions [SOI, 



Vpi,{ql,iujs) 



(8) 



where q\ = 2k'j{l — cos cj)), is the momentum exchange 
between the interacting particles along the Fermi sur- 
face. At r = 0, assuming momentum transfer occur 
around Fermi momentum kf, we can expand Eq. ([S]) as 

Vphikl,0) = - ., -1,0,1,. and the di- 

mensionless effective interaction between the fermions in 
angular momentum channel L is given by 



Al(0,0) = Gm 



2tt 



exp[iL(j)\d(l)/2TT 



FERMIONIC SELF-ENERGY INCLUDING 
VERTEX CORRECTION 



In this section, we consider the fermion self-energy 
in the normal state due to the interaction between the 
fermions and dressed Bugoliubov quasi-particles. In do- 
ing so, we explicitly take into account the momentum de- 
pendance as well as the retardation of the phonon prop- 
agator. More importantly, we also go beyond Migdal's 
adiabatic limit and take into account the effect of ver- 
tex correction. This kind of non-adiabatic correction is 
introduced to electron-phonon system in metals in Refs 
[sH - fssj in the limit of T — 0. Our main result is that as 
temperature gets lower, the vertex-corrected mass renor- 
malization function Z{T) also gets smaller. 

The unperturbed Green's function for the fermions is 
given by Go{k±, iujn) = l/{iujn — £,{k±)), where the Mat- 
subara frequency uj„ = {2n + l)TrT, n being integer. The 
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COS 
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(9) 
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As we are considering single component fermions, pairing 
will occur in the odd angular momentum channels. 



FIG. 4: Figure of fermion self-energy in normal state including 
the first order vertex correction. The solid hue denotes the 
fermion propagator while the thick dashed line denotes the 
dressed phonon propagator of Eq. ((6)1 




FIG. 3: Figure of effective fermion-fermion interaction 
strength A_t(0,0) in the angular momentum channels L — 
0(dashed line), L = l(black line),= 3(red line),= 5(blue 
line) as a function of the fermion dimensionality parameter 
rj. We fixed mf/mb = 53/52, gsd = 4, Sg/STrgdd = 0.6, and 
Qm = 0.5. 



normal self-energy T,n{k±), including the first order ver- 
tex correction, as shown in Fig. |4l is given by 



E„(fc_L,ia;„) = ^^^(A:^, ia;„) -f I]^(fc_L, iw„), (10) 



where I]^j(/c^, ia;„) comes from the first diagram in Fig.|4] 
and Y^'!^{k±,iLtJn) arises from the second diagram which 
denotes the vertex correction. Also, 



5 



m,q1 

T:l{k±,iuJn) = ^ \j{kj_ - ql_)\'^\j{k± - pl)\'^D{uJi - ujn,k± - pl)D{uJm - U!n,kj_ - ql) (11) 

m,l,q1,p1 

X Go(ia;m,<7l)Go(ia;;,pl)Go(i(w; - w„ - fc^ 

(12) 

I 

By taking the average over Fermi energy(|fcj_| « \q±\ w and the angles {k±,qj_) — (j)' , {k_L,'p\_) = 6, ^ — (j)' — 

kf), and denoting ^{k±) = h^k\l2m~ e/ x, £_{pl_) = ^- Using Euler-Mclauren summation formula, we can 

h^kl/2m - e/ = x', and ^(pl - k^ql) ~ = x' + 2efa, transform the sum over I in Eq. (HH) to integral as, 
where 

a — 1 — cos((/)') + cos(6') — 003(7) (13) 



P{x',m,T,(l)',0) = - D{uji - ujn, fc_L - pl) 



(14) 



X G{iujrn, ql)G{iuji, pl)G{i{uji - a;„ +a;„i),pl - fc_L +gl) 

duj 1 1 

J 27r + fl'^ iuj — x' + inT {uj + 2rmrT) — x' — a + ittT 



in 



2(2m7rT - ia) 



e{x') 



i^x') 



— ittT VL — X 



[x + a) 



iirT 

—x' — a) 



fl — 2irmrT + x' — inT + a fl + 2irm:T — x' + iirT — a 



(15) 



here 0{..) is the Step function. We can rewrite E„(za;n) 
as 

E„(iw„) = x{i^n) + iujn[l - Z{iujn)], (16) 

where x(iw„) is the real part of I]„(iti;„) and Z{iuJn) 
is known as the mass renormalization function. x(7rT) 
shifts the bare fermion dispersion energy ^{k±). As 
we are considering strong-coupling superfluidity, we only 
consider the term corresponds to n = 41]. Then the 
mass renormalization function Z{ttT) and x(7rr) in the 



conventional Eliashberg form is given by, 

ZinT) = l + XziT), (17) 
xiirT) = ^sign(c^„)Im[P4m,r)], (18) 

m 

with an effective coupling constant 

Az(r) = Ao(0,0)-P,(r), (19) 
and P,{T) = ^„ sign(c^,„)Re[P,(m,T)], where 



I 

Xl{s,T) = - / \-fi^)fD{ujs,<P)cosL(j)dcj)/2n, 
Jo 

Pv{m,T) = r dx' r ^d0\^{cl>')\'\j{e)\'P{x',m,T,<p\e), 

(20) 



and Re[/], Im[/] denotes the real and imaginary part of the function /. In Fig. [5] wc plot Pv{T) as a func- 
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tion of temperature for various values of dimensionality 
ij. From the expression of Xl{s,T) in Eq. ([20]) . we no- 
tice that Al(0, 0) is always positive as shown in Fig. |3l 
Py{T) / Xo{0, 0) is the expansion parameter for the pertur- 
bative scheme. The correction Pv{T) is always found to 
be positive, which reduces the coupling strength Xz{T). 
With decreasing temperature, Pi,(T)/Ao(0, 0) increases 
before saturating. This saturation value increases with 
decreasing dimensionality as seen in Fig.O For smaller 77, 
P„(T) Ao(0, 0) becomes higher than one for lower temper- 
ature, invalidating any perturbative calculation in that 
low temperature region. 




the same parameters as Fig. [5j We find that x(T) <C 
Ao(0, 0) even in the temperature region with high Py(T). 
Henceforth, we can neglect the effect of energy shift of 
the fermions due to the smallness of x{T). 

VI. SELF-ENERGY IN COOPER-PAIR 
CHANNEL AND TRANSITION TEMPERATURE 

In this section we look at the fermionic self-energy in 
the Cooper-pair channel, Yjs(k±,iu}n) as represented di- 
agramatically in Fig. [71 close to transition temperature 
Tc- Then we find transition temperature within strong- 
coupling limit considering the terms n = 0,-1 [33 |. The 
main results of this section are : i) Depending on the an- 
gular momentum channel, the vertex-corrected interac- 
tion strength can increase as a function of dimensionality 
and temperature, ii) The solution of modified Eliashberg 
equation supports Tc ~ .lef in the strong-coupling limit 
for p—,f—, and h — wave order parameters. 



FIG. 5: Plot of the vertex correction in as a function of 
for different values of the dimensionality parameter 77 = 
1.5(the dotted line), 77 — l(the solid line), rj — .6(the dashed 
line). We fixed mf/nib — 53/52, = -6,5 — .5, and 

gsd = 3.8. 



Next, we calculate the real part of the fermion self 
energy xi^T). Figure [6] shows a such generic case for 




FIG. 6: Figure of x(r) as a function of ^ for different values 
of the dimensionality parameter -q — 1.5(the dotted line), — 
l(the solid line), rj = .6(the dashed line). We fixed mf/mi, = 



53/52, 



3g 
4TSdd 



— ■G,Q = .5, and gsd — 3.8. 



+ I 



+ ^ 



FIG. 7: Schematic diagram of fermion self-energy in the 
Cooper-pair channel channel including the first order vertex 
correction. The solid line denotes the fermion propagator 
while the thick dashed line denotes the dressed phonon prop- 
agator of Eq. ((6]). The first diagram denotes the direct inter- 
action, while the second and third diagram denote the first 
order vertex corrections, and the third diagram denotes the 
cross-term. 



For T > Tc, the fermion self-energy in the Cooper-pair 
channel, I]s(/c^), is given by 

Y.s{k±,iujn) = T,l{k±,iujn) + S^(fcj_, ii:j„) -|- iw„), 

(21) 

where (fcj_, za;„) comes from the first diagram in Fig. [71 
the second and third diagram gives equal contribution 
and denoted by S^(/cj^, iaj„) while the last diagram de- 
noting cross interaction is given by S]^(fcj^, iujn), and 



S^(fcj.,iw„) = -T ^ |7(fc_L - ql)|^i:)(a;„i - w„, fc_L - gl)G'(ia;m, gl)G(-ia;m, -i_L)Ss(9l,«'^m), (22) 
S^(fcx,iw„) = 2r^ ^ \j{kj_-ql)\'^\j(kj_-pl)\'^D{uj,n~oJri.,kj_~ql)D{uJn-(^i,k±-ql) 

X G{iuJi,pl)G{i{uJi -ojn +ujjn),pl - k±_ + ql_)G{iuj„i,ql)G{-iujjn, -q±)'^siql,ii^rn), (23) 
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X G{iuji,pl)G{i{uJi - w„ - oj,n),pl 



We take the average over Fermi energyd A; _L I « \q±\ « fc/), 
and denote ^(pl — — q±) — ^ = x' + '2ef(3, where 

(3 — 1 + cos{4>') — cos(0) — cos(7). 

We define the superfluid order parameter in the usual 
way, A(«a;„) — 'Z,a{iuJn) / Z{iuin)- Due to single- 
component nature of the fcrmions in the mixture, the 
order parameter can be expanded in odd partial wave 
A(ia;„,(/>) = X]l=...,-i,i,... ^i(*'^n) expiii^]. In the 
strong coupling limit we are interested in the terms n = 
and n = — l[4l|. Assuming the order parameter to be 
even function of frequency, A(7rT) = A(— ttT), we get 



il)pD(a;„ - uji,k±_ -pl)D{uJrn ~ oJi-.ql - pl) 

(24) 



from Eq. (Ell), (E2I), 

Ai(7rT)Z(7rT) - A^(0, T)Ai(7rr)+Af (-1, T)Ai(-7rT), 

(25) 

where 

\t{m,T) - AL(m,T)-2P,(m,T)-P,(TO,r), 

(26) 

XLim^T) originate from Yi\{iujn), whereas P^{m,T) 
comes from vertex corrected self-energy I]g(ia;„). 
Pc{m^T) results from the contribution the cross-term 
Ss(iw„), 



P^{m,T) = 2TIm 



E l7Wl'l7(^ - - ^')D{m - I, 0') 



tan 



{21+1)t:T 



— tan 



i/3 (2/ - 2m - l)7rr 



;3-F2i(m-h l)7rT_ 



(27) 



Eq. (pS)) is similar to the Eliashberg equation [41] in the 
strong coupling limit with additional vertex-corrected in- 
teraction strengths. Apart from the dependance on tem- 
perature, the effective interaction strength in the angular 
momentum channel L, (T) , is also a function of dimen- 
sionality r] and the boson- fermion mass ratio rrif/mb. As 
noticed in Eq. (|20|) . Pv{0, T) is positive, thus reducing the 
effective interaction strength in the Cooper-pair channel. 
But the correction arising from the cross term Pc{m,T) 
can be negative or positive, as shown in Fig. \8\ Also the 
sign of Pc{m,T) depends on the particular angular mo- 
mentum channel under consideration. But in general, for 
smaller 77, Pc{m,T) becomes positive, in turn reducing 
the interaction strength, X^{m,T), in the Cooper-pair 
channel. But higher value of 77 changes Pc{m,T) to neg- 
ative, which enhances X'^{m,T). We also find out that 
with increasing fermionic mass, Pc{m,T) becomes nega- 
tive for lower values of 77, thus enhancing the interaction 
strengths in the various angular momentum channel. The 
reason behind this is that Pc{m,T) depends on the ra- 
tio between the position of the roton minimum fco ^nd 
Fermi momentum k()/2kf. With high fermionic mass, kf 
increases, resulting in lower ratio fco/2fcy. This in turn 
makes the cross term more negative. From this we infer 




0.5 1 1.5 2 

1 



FIG. 8: Plot of cross-term interaction Pc(0, e/) as a function 
of 77 for mrif /vrih — 53/52, ttT — fi, gsd ~ 4.0, Gbf = 0.5 and 
= 0.6. The black, red and blue lines correspond to p, 
f, and /i-wave angular momentum channel respectively. 



that it is better to use fermions with higher mass to get 
a higher interaction strength A^(to,T). 

Next we look into the dependance of A^(0,T) on di- 
mensionality 77. First we fix ttT — £f and plot for var- 
ious value of 77 n Fig. [HI We see that the magnitude of 
A^(e/) in different angular momentum channel changes 
as 77 varies. Also different angular momentum channel be- 
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FIG. 9: Figure of effective interaction X^{tvT = e/) as a func- 
tion of T] for rrif/mb = 53/52, ttT = e/, gsd — 4.0, Qht ~ -5 
and — = .6. The black, red and blue lines correspond to 

Af (r),!^ (T), (T) respectively. 



comes dominant depending on the dimensionality. This 
qualitatively resembles the situation in section 3, where 
without the vertex correction, we find that different in- 
teractions becomes dominant for different dimensionality. 

Next we find the transition temperature for fermionic 
superfluidity within perturbative scheme as long as < 
T* where T* is the temperature for which the pertur- 
bative scheme becomes invalid. In the situation of more 
than one solution of Eq. (|25|) . we have taken the tran- 
sition temperature to be the one corresponding to the 
highest temperature. In the case of Eq. (P5|) having no 
solution for Tc, vertex-corrected strong coupling super- 
fluidity is not possible and to find the transition temper- 
ature we need to consider the full Eliashberg equation 
for all values of n. This regime is not considered 
in this paper as we are only interested in the high tem- 
perature limit of the transition temperature. In Fig. [TU] 
(a) and (b) we plot the solution of the Eqs. ([TT]), ([25l) . 
In the case of obtaining multiple solutions for transition 
temperature corresponding to different angular momen- 
tum channel, we considered the channel with maximum 
transition temperature to be the solution. In general we 
find that with lowering dimensionality, the nature of the 
superfluidity changes from p, to h to /-wave. Also we 
obtained transition temperature Tc in the order of Aej 
or more in each angular momentum channel. By compar- 
ing Fig. [To] (a) and (b) we notice that, for higher fermion 
mass, we flnd that much lower value of rj can be attain- 
able without the perturbation becoming invalid than for 
a lower fermionic mass. 



VII. CHIRALITY AND NON-ABELIAN 
ANYONS 



0.8 1.2 



(a) 




(b) 

FIG. 10: Critical temperature Tc, in the units of Fermi energy 
6/ , as a function of rj. The black, red and blue lines correspond 
to transition temperature of p-wave, h-wave, and /-wave or- 
der parameters respectively with parameters . a) For param- 
eters nif/mt — 1.7, and dipole strength g^^ ~ 2.5. b) For 
parameters mf/rrib = 53/52, and dipole strength g^d = 4.0. 



are bounded and non-Abelian in nature. By consider- 
ing the superfluid gap equation at low temperature, the 
gap is maximum when the order parameter breaks time 
reversal symmetry. From now on we assume that the or- 



1 L 



iL0 



der parameters are denoted by = Ao(r) 

where = kcosO, ky = ksmO and Ao(r) is the cen- 
ter of mass amplitude of the Cooper pairs, with f being 
the center of mass coordinate of the pair. For vortex 
state Ao(r^ cab be approximated as: i) Ao(r) — Q,r<^ 
and ii) Ao{r) ~ Aq exp{i(j)), r > ^, where r = yaJ^-l-y^, 
tanc/) — y/x. ^ is the size of the core of the vortex. 
The vortex state of the p-wave superfiuids always has a 
zero-energy bound quasi-particle state IH, [s^ . Now 
we discuss the asymptotic solutions for the zero-energy 
bound state for /- and h-wave order parameters. The 
quasi-particle states in a single vortex can be found in 
the limit of large distance from the vortex core by solv- 
ing the Bugoliubov-DeGennes equation. 



In this section we briefly discuss the properties of 
quasi-particle excitations inside a vortex for different in- 
ternal symmetries of the order parameter. By solving 
Bugoliubov-DeGennes equation in the limit of large dis- 
tance from the core of the vortex, we show that the zero- 
energy solutions in the chiral p—, /— , h— wave states 



HqUl + (-«)"Tre 



V 



■\L^ ~i(t>/2 

k^ 



20 



d id 
dr r dcf) 



d i d 
dr r d(j) 



Eur 



-ii>l2 
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(28) 

where E is the energy of the quasi-particles with ampli- 
tudes ul,vl- We particularly look for zero energy solu- 
tions with bounded ul,vl with the property ul — v]^ 
[ssj . For r — >■ oo, we can neglect the terms scaled as 
in Eq. Then the solution of Eq. ([28]) with different 

orbital symmetries read, 



exp 
cxp 



cxp 



6m f A I 



■r 

1/3 



(29) 



The zero-energy solution for each odd-wave parameter 
corresponds to different angular momentum channel of 
the quasi-particles inside a vortex core. These results 
can also be carried out by applying the "index theo- 
rem" 40] . For temperature smaller than the energy gap 
Ao/e/, only the zero energy mode is occupied. The 
quasi-particle operator in that situation is written as 
7l = f (Pr(uL(r )c' (r) + VL(r)c(r)), which acts as Majo- 
rana fermion |6l. Il2l l38l|. 71, obeys non-Abelian statistics 
and can be used for quantum computing 13]. In or- 
der to perform quantum computational tasks, existence 
of several well separated vortices is necessary. We can 
assume in the weak coupling limit, Ag = /3e/, where 
/? is a constant usually less than one. Then substitut- 
ing Aq in Eq. (P^ . we get ui oc exp(— /3fc^r/2), U3 c>c 
exp(— /c/r/3/3), U5 cx exp(— fc/r//?-'^/'^). For r '3> for 



smaller /?, 1*5 has smaller tail than and ui. Thus 
non-overlapped states can be achieved more easily in su- 
perfluids in i = 5 and L ~ 3 channel than for p-wave 
channel. 



VIII. CONCLUSION 

Summarizing, we studied boson-induced superfluidity 
of fermions in a mixture of dipolar bosons and single- 
component fermions. A system is proposed where con- 
ventional pairing mechanism gives rise to different exotic 
internal structures of the Cooper pairs with strong in- 
teraction in respective angular momentum channels. We 
find that vertex-corrections play an important role for 
superfluidity in this mixture and results in high value 
of transition temperatures. We like to stress that the 
high transition temperatures are a result of the inclu- 
sion of vertex correction and cross interaction within the 
Cooper-pair channel. Importantly, we find that by de- 
creasing ?7, we can generate exotic superfluids with p, 
f and h-wstve internal structure. Excitations in these 
types of superfluids breaks time-reversal symmetry and 
supports quasi-particles with non-Abelian statistics. 
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